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O [ Abstract 
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■* — . The most precise test of the post-Newtonian 7 parameter in the solar system has 

qh| been achieved in measurement of the frequency shift of radio waves to and from the 

Cassini spacecraft as they passed near the Sun. The test rehes upon the JPL model 
bX)' of radiowave propagation that includes, but does not explicitly parametrize, the 

impact of the non-stationary component of the gravitational field of the Sun, gen- 
ts^ I erated by its barycentric orbital motion, on the Shapiro delay. This non-stationary 
j_j ' gravitational field of the Sun is associated with the Lorentz transformation of the 
^ metric tensor and the affine connection from the heliocentric to the barycentric 
frame of the solar system and can be treated as gravimagnetic field. The gravimag- 
netic field perturbs the propagation of a radio wave and contributes to its frequency 
shift at the level up to 4 x 10~^^ that may affect the precise measurement of the 
parameter 7 in the Cassini experiment to about one part in 10,000. Our analysis 
suggests that the translational gravimagnetic field of the Sun can be extracted from 
the Cassini data, and its effect is separable from the space curvature characterized 
by the parameter 7. 
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1 Introduction 



According to general relativity and other alternative theories of gravity, pho- 
tons are deflected and delayed by the curvature of space-time produced by 
any mass distribution. In the gravitational fleld of a static spherical mass 
the bending and delay are proportional to 7 -|- 1, where 7 parametrizes the 
curvature of space and measures the degree of deviation of a gravity theory 
from general relativity [1] . The most precise measurement of 7 has been re- 
cently achieved in the Cassini experiment [2] by fltting the Doppler data to the 
general- relativistic model used in the JPL Orbit Determination Program [5]. 
The ODP calculates the time delay in the heliocentric frame and transforms 
it to the barycentric frame of the solar system by using spatial translations of 
participating bodies from the heliocentric to the barycentric frame. 

We demonstrate in this Letter that the spatial translations used in the ODP are 
equivalent to performing a Lorentz transformation of the solar gravitational 
fleld from the heliocentric (H) to the barycentric (B) frame. In the heliocentric 
frame photon moves in a static gravitational fleld of the Sun that is fully 
deflned by the metric components Qqq and g^j , while g^^ = 0. The Lorentz 
transformation of the gravitational fleld to the barycentric frame brings about 
the gravimagnetic, (7^ 7^ 0, component of the metric tensor [1]. Moreover, 
it makes all components of the barycentric metric time-dependent, so that 
its partial time derivatives do not vanish. This non-stationarity of the solar 
gravity in the barycentric frame is associated with the gravimagnetic fleld 
caused by the translational current of Sun's matter [5] lJ- 

This Letter shows that the barycentric orbital velocity of the Sun gives rise to 
gravimagnetic shift of frequency, zqm ~ C(qVq/c, where a© is the solar gravi- 
tational deflection of light, Vq ~ 15 m/s is the barycentric velocity of the Sun 
during the time of the Cassini experiment 2002, and c is the speed of light. For 
the light grazing the solar limb a© = 8.5 x 10~^ rad, and zqu < 4.3 x 10~^^. 
Therefore, detection of the gravimagnetic frequency shift zqm is to be possible 
in ranging experiments like the Cassini having precision 10^^^, or better. Dur- 
ing the time radio signal travels from the spacecraft to the Earth the motion 
of the Sun can be approximated as uniform along a straight line as the solar 
acceleration, a©, produces relativistic effect that is negligible for the Cassini 
but may be detected in future experiments [7]. This acceleration-induced fre- 
quency shift Zacc. ~ aQaQR/c^ ~ aQ^R/ Rq){vq/cY <ti lO^^'', where R is 
the distance spacecraft-observer, Rq is the solar radius, and we have used 
Qq ~ Vq/Rq as the radius of the solar orbit around the barycentre is about 



^ Notice that GP-B experiment [6] measures the gravimagnetic field caused by 
Earth's rotation mass-current that can not be generated by the Lorentz transfor- 
mation. 



Rq. 

In what follows, the Greek indices a, (3, ... run from to 3, the Roman indices 
i,j, ... run from 1 to 3, repeated indices mean Einstein's summation, and bold 
letters A = {A^),B = (5*), etc. denote spatial vectors. A dot between two 
spatial vectors, A ■ B, means the Euclidean scalar product, and the cross 
between them, A x B, means the Euclidean vector product. The absolute 
value of a vector A is denoted A = \A\. 



2 Gravitational Shift of Frequency in the Heliocentric Frame 



The heliocentric frame is denoted by X° = {cT,X), where T is time and 
X is the space coordinates with the Sun located at the origin. Spacecraft 
emits a radio signal at time Tq at point Xq, and observer (tracking station) 
receives this signal at time Ti at point Xi. The radio signal propagates along 
the direction defined by a unit vector K = R/R, where R = Xi — Xq, and 
R = \Xi — Xo\. The total time of propagation of the radio signal from the 
spacecraft to the observer (with plasma dispersion eliminated) is 

ri-To = - + AT, (1) 

c 

where AT is the Shapiro gravitational time delay [8] 

GMq, fRo + Ri + R 



G is the universal gravitational constant, Mq is mass of the Sun, Rq = \Xq\, 
Ri = \Xi\. Equation ([2]) is valid in the post-Newtonian approximation and 
includes all relativistic terms of the order of 1/c^ in addition to the first term 
in the right side of equation ([T]). In the precise experiments like the Cassini 
some terms of the higher order of 1/c^ must be included to equation ([2]). They 
were derived, for example, by [9] and [10]. The most important of these higher- 
order terms are included to the ODP but are omitted here because we focus 
on the discussion of the velocity-dependent, gravimagnetic effects while the 
terms of the higher order are static. 

The ODP computes Doppler shift by differencing numerically the time delay 
equations ([I]), ([2]) at the beginning, Tq, and the end, Ti, of the integration 
time. This procedure needs better theoretical understanding and is described 
here analytically in the approximation of a small impact parameter D of the 
radio link that is sufficient for our subsequent discussion. In this approximation 



equation ([2]) is simplified 



AT^,1H-.)^,„(1^), (3) 



where Li = (K ■ Xi) = ^Rf - D^, Lq = -{K ■ Xq) = ^ Rl - D^, D = 
K X (Xi X K) = K X (Xq X K), D = \D\, and we notice that Li + Lq = R. 
In the hehocentric frame, only the relative motion of the spacecraft and the 
Earth produce a spatial change in the impact parameter leading to a temporal 
change in the delay AT which re-calculates to a change in observed frequency 
z/i of the radio signal emitted from the spacecraft with frequency z/q. The 
frequency change is z = z/i/i^o ~ 1 = zq + zi + Zg,., where Zq and Zi are 
corrections accounting for transformation from time T to the proper time r 
at the spacecraft and observer positions respectively, and Zgj- = d{AT)/dTi is 
the delay due to the variation in the strength of the gravitational field along 
the radio signal path due to the motion of observer, satellite, and/or the 
gravitating mass [TT] . Only Zg,- is important for analysis of the gravimagnetic 
effect of the barycentric motion of the Sun in precise determination of 7. Its 
calculation yields 






V = -7^(Vi ■ an) - ::^(Vo ■ «^) + Sz^, . (4) 



Here Vi = dXi/dTi and Vq = dX^/dT^ are the heliocentric velocities of 
observer and emitter with respect to the Sun, residual terms 



^ 2c 



^{K.V,)-^{K.V,) 

Ml 1^0 



^' K-{Vi- Vo) {Vo ■ cxh) , (5) 



c^ 



where an = \(^h\, 

"/f = «0^— -^-D (6) 

is a displacement vector describing magnitude and direction of the gravita- 
tional bending of the radio wave's trajectory, D = D/D. The asymmetry 
between subscripts 1 and in the last term of equation ([5]) is because of the 
derivative dTo/dTi = l — c~^K-(Vi — Vq) that appears any time when one dif- 
ferentiates functions depending on Xq = X(To(Ti)) as dX^/dTi = VodTo/dTi. 



The first two terms in the right side of equation (j4]) were obtained by [12]. How- 
ever, the analytic theory of the Cassini experiment should also incorporate the 
term Sz^j. given in equation (^. Indeed, the r.m.s. value of data points in the 
Cassini measurement corresponded to a one-way range rate of 6V = 2.2 x 10~^ 
m/s [2]. Therefore, all terms having amplitude 6V/c ~ 10~^^ -7- 10^^^ must be 
kept in the analytic approach. The ODP calculates Sz^j. numerically [3]. The 
analytic treatment of the gravitational Doppler shift provides deeper theo- 
retical insight to the ultra-precise Doppler measurements in space navigation 



that is used for prediction more subtle relativistic effects [TT]|13fl4|[T5] which 
are feasible for detection in future missions [7]. 



3 Gravitational Shift of Frequency in the Barycentric PVame 

3.1 The Lorentz Transformation 



The barycentric frame is denoted by x" = (ct, x) where t is a coordinate time 
and X are space coordinates with the origin located at the barycentre of the 
solar system [TB]. Let x^it) denote the spatial coordinates of the Sun at time 
t and its orbital velocity '^©(t) = dx^/dt. Radio signal is emitted at Xq at 
time to; and received by observer at Xi at time ti. Because we neglect the 
solar acceleration, a©, its coordinate can be approximated by uniform and 
rectilinear motion 

XQ{t) =XQ + VQ{t- tA) + O [aQ(t - IaY] , (7) 

where Xq = XQitA) and Vq = VQ{tA) are respectively the barycentric coordi- 
nate and velocity of the Sun taken at a fiducial instant of time, tA, which is 
in the range between to and ti. 

In what follows we assume that relativistic transformation of the solar gravi- 
tational field from the heliocentric to the barycentric frame can be performed 
with the matrix of the Lorentz transformation A^ having a standard special- 
relativistic form [T7] 



A°o = r , A° = A*o = -r/3' , A^ = 6'^ + (L - 1) ^'^' 



P' 



where T = v^l — (3'^, and /3* = Vq/c <^ 1. This matrix transforms coordinates 
as follows 

X" = A^ {x^ - 4) , (9) 

where Xq = (ctA, Xq). The Lorentz transformation transfrorms the metric ten- 
sor and the afiine connection, and make them time-dependent in the barycen- 
tric frame of reference. All time-dependent terms that are proportional to the 
barycentric velocity of the Sun, Vq, are called gravimagnetic [^j and are as- 
sociated with the gravimagnetic effect in the Doppler shift. Precise equations 
describing the Lorentz transformation of the gravity field and their connection 
to the propagation of light are given in [1]. 

As we are interested only in the gravimagnetic shift of frequency of the first 
order, it suffices to approximate equation (^ by making use of the Taylor 



expansion of A^ with respect to f3 and to work with a few first terms. This 
reduces equation ([9]) to the post-Gahlean transformation 



X = x-XQ-VQ{t-tA) + 0{vl/c'), (10) 

T = t + 0{v^/c'). (11) 

It is important to emphasize that the term proportional to Vq in equation flTUl) 
is, effectively, a linear order with respect to Vq/c, as it is obtained from the 
contraction of A*q from equation ([8]) with time coordinate x^ = ct. Therefore, 
the Lorentz invariance of the gravity field reveals itself already at the level of 
the post-Galilean transformation (fTOl) that can be tested experimentally with 
the Cassini data. 



3.2 Transformation of the Gravitational Shift of Frequency 



Transformation of the frequency shift equation (jlj) to the barycentric frame 
requires transformation of the heliocentric distances Rq, Ri, R, the impact pa- 
rameter D, and velocities Vq, Vi to their barycentric counterparts. By making 
use of equations flTUl) . fITT]) we obtain 

Xi = ri=a;i-a;0(ti) , (12) 

Xo = ro = xo-XQ{to), (13) 

V, = v,-Vq, (14) 

Vo = t;0-t7Q, (15) 

where XQ(to) and XQifi) are defined by equation (jTj), vi = [dxi{t) / dt]t=t-^ and 
Vq = [dxQ{t)/dt]t=tg are the barycentric velocities of observer and emitter of 
the radio signal respectively. 

Equations (fT2l) . (TT3l) are the spatial translations used in the ODP for calcu- 
lating the barycentric time delay [3] . They allow us to derive transformations 
of the distance R and the unit vector K. Specifically, 



R=r-'^r, (16) 

c 

K = k-^, (17) 

c 

where r = Xi — Xq, r = \r\, the unit vector k = (xi — a3o)/|a;i — Xq\ defines 
direction of propagation of the radio signal in the barycentric frame, and v^ = 
k X {vq X k) is the Sun's velocity projected on the plane being orthogonal to 



k. Transformation of the heliocentric vector D is given by equation [18 



D = d+{d-VQ)k + {k-ri)v^, (18) 

where d = k x (ri x k) is the barycentric impact parameter vector of the 
hght ray taken at instant ti. Using equation (TT5]) the deflection angle (Xh is 
transformed as follows |18] 



cxh = olb + {ctB ■ VQ)k + 



{k ■ rO 



asVQ -2(d- Vq)cxb 



l+7i?o^ 



d 



(19) 
(20) 



where as = |ckb|, and the unit vector d = d/\d\ 



Substitution of the hehocentric quantities (fT2|) - (pO|) to the frequency shift 
equation (jl]) yields 



^0 / N 

— {vi -olb) 
cr 

observer shift zq 



— (vo ■ olb) + -(*^0 ■ olb) 

cr c 



+5zl 



satellite shift zg 



gravimagnetic shift zqui 



(21) 

where Iq ^ —{k ■ r^), li = {k ■ ri), and 6z^^ denotes the residual terms with 
magnitude up to 10~^^ [H]. In the Cassini experiment the distances /q — 7.4 
AU, h ^ 1-0 AU, and r = 8.43 AU. 



3.3 Parametrization of the Shapiro Delay in the Barycentric Frame 



Linearized transformations ( IT2l) and ( IT3l) of the Shapiro time delay are used in 
the ODP |3] but they are not parametrized to track down the gravimagnetic 
effect of the solar barycentric velocity explicitly, which may be also important 
for other, ultra-precise solar system experiments [7]. The most simple way to 
introduce such a parametrization is to label velocity of the Lorentz transfor- 
mation by replacing /?* —>■ e/3\ where e is a new fitting parameter [1]. It is not 
equivalent to the PPN parameter ai as ai is associated with the "preferred- 
frame" effects [1] which are absent in general relativity, while e marks genuine 
gravimagnetic effects induced by the Lorentz transformation of gravitational 
field [H]. The Cassini data are accurate enough to measure e, thus, providing 
a test of the existence of the gravimagnetic field caused by the barycentric mo- 
tion of the Sun. Gaussian distribution of the resiuduals of the post-fit Cassini 
data [2] signals that e is close to unity. However, it is tantalizing to obtain its 
precise statistical estimate that will provide a new test of general relativity at 
higher technological level. 



The suggested e-parameterization recasts equations (fT2l) . (fT3l) to the following 
form 



Xi = ri=xi- XQ{eti) , 
Xo = ro = Xq - XQ{eto) , 

so that equation (|2T|) becomes 



'gr 



^0 / N 

cr 



cr 



-{vq -olb) + 6z (e) 
c '^ 



(22) 
(23) 



(24) 



making evident that e traces that part of the Doppler shift which is solely due 
to the barycentric orbital motion of the Sun. The partial derivatives of the 
Shapiro delay with respect to parameters 7 and e are 



dAT 

dAT 
~d7~ 



Je=l 



7=1 




k 



(ri + ro)2 - r 



(25) 



(ri + r2) ,(26) 



where fci = Vi/ri and fco = '^o/^o- Statistical meaning of partial derivatives 
fl25l) and fl26l) is that they are components of the covariance matrix of the fitting 
parameters. The fact, that they have different functional dependence, points 
out that parameters 7 and e are uncorrelated and, hence, can be separated 
with confidence in the data fitting procedure. 



4 Discussion 



The Cassini experiment was conducted in June 2002 when radio waves from/to 
the spacecraft had passed at the minimal distance d = 1.6 Rq from the Sun. 
The peak value of [-Zgrjpeafc — (^Q{RQ/d){v^/c) ~ 6 x 10""'^° is caused by the 
orbital velocity v^ of the Earth 
about five orders of magnitude smaller [2] . 



while the measurement errors of Zg^ were 



We have used JPL ephemerides to calculate the barycentric position and ve- 
locity of the Earth and the Sun at the time of the Cassini-solar conjunction. 
One of the consecutive configuration of these bodies is shown in Fig. [TJ It 
varies rapidly as all three bodies are moving with respect to each other, thus, 
leading to the change of the impact parameter of the radio link. The solar 
barycentric velocity is about 14.4 m/s, and the angle between the direction 
of Vq and the vector d is about 16°. This angle changes its value to ^ 164° 



after the satellite passes on the other side of Sun's limb, thus, changing the 
sign of the gravimagnetic effect. 

These values were used for comparative estimation of the gravimagnetic term 
in equation fl2T]) . The Doppler residuals in the experiment run by [2], were 
obtained from 9-parameter fit on the calibrated data with the solar corona 
and the interplanetary medium dispersions removed with double-frequency 
observations. The r.m.s. value of the N = 1094 data points (relative to an 
8.4 GHz carrier) is 1.2 x 10~^ Hz, corresponding to a one-way range rate of 

2.2 X 10~^ m s~^ or for a relative frequency shift 6z = 0.7 x 10~^^. It means that 
the average error of a single measurement is az = y2N5z ^ 3.2 x 10"^^. The 
gravimagnetic frequency shift caused by the barycentric orbital motion of the 
Sun amounts to the peak value of [^civilpeafc = ij /lo)iyQ/v^) cos (t)[zg-c]peak — 

3.3 X 10~^^ that is comparable with the r.m.s. of a single measurement but 
a comprehensive number of measurements is required to extract zqm out of 
noise. In order to estimate the number of the measurements being available 
for this goal from the existing data, we notice that the gravimagnetic term 
^GM ~ (%/c)aB reaches the Doppler data r.m.s. value for d ^ 30 Rq that 
corresponds to 20 days of observations around the solar conjunction giving 
the number of data points n ^ 800. This should be statistically sufficient to 
measure the gravimagnetic effect caused by the barycentric motion of the Sun. 

Equation (1211) suggests that the gravimagnetic term -Zgm, if it were ignored 
for whatever reason, would affect the measured value of the parameter 7. The 
bias in 7 is approximately given by the ratio [zGyi\peak/[zgT\peak < 5 x 10~^ 
which should be further relaxed by a factor of 3 -^ 5 for 2;gm is extracted 
from a smaller subset of the entire data [fj. This value is larger than 7 — 1 = 
(-1-2.1 ±2.3) X 10~^ given by |2] pointing out that the gravimagnetic field due 
to the solar barycentric velocity Vq can be measured though its extraction 
may represent difficulties (Bertotti 2006, private communication). 

This Letter suggests to disentangle the measurement of the PPN parameter 
7 in the Cassini experiment [21121] from the effect of the gravimagnetic field 
caused by the barycentric orbital motion of the Sun. The orbital motion must 
be properly parametrized in the delay equation ([2]) by making use of the 
translations fl22l) . fl23|) where e is a fitting parameter. This paper calculates 
the gravimagnetic frequency shift in one-way radio link. In real experiment 
two - up and down - radio links are used. Analytic equations of the Doppler 
shift for the two radio links were derived in [TT] and they substantiate our 
conclusions. 

The translationally-invariant gravimagnetic field can be measured more ac- 
curately in space missions LATOR [22] and ASTROD [23]. The space inter- 



have arrived to similar numerical estimate by analyzing solar-system experi- 
ments in the framework of the Standard Model Extensions. 



ferometer SIM [21] and ground-based Square Kilometer Array [25] can also 
achieve this goal [26]. 
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Fig. 1. Relative configuration of the Cassini, the Sun and the Earth by the time of 
solar conjunction on June 21, 2002 in the barycentric reference frame with origin 
located at the centre of mass of the solar system. The barycentric velocity of the 
Sun is shown by (red) arrow that makes an angle (p ~ 16° with respect to the im- 
pact parameter of the radio-link connecting the Cassini spacecraft and the Doppler 
tracking station on the Earth. This angle changes to (/> ~ 164° after the radio-link 
gets on the other side from the Sun. 
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